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Screening effects of eiectrons on inhomogeneous nuciear matter, which inchides spherical, slablike, 
and rodlike nuclei as well as spherical and rodlike nuclear bubbles, are investigated in view of 
possible application to cold neutron star matter and supernova matter at subnuclear densities. 
Using a compressible liquid-drop model incorporating uncertainties in the surface tension, we find 
that the energy change due to the screening effects broadens the density region in which bubbles 
and nonspherical nuclei appear in the phase diagram delineating the energetically favorable shape 
of inhomogeneous nuclear matter. This conclusion is considered to be general since it stems from a 
model-independent feature that the electron screening acts to decrease the density at which spherical 
nuclei become unstable against fission and to increase the density at which uniform matter becomes 
unstable against proton clustering. 



PACS numbers: 97.60.Jd, 26.60.+C, 64.75. +g, 21.65.+f 



I. INTRODUCTION 

At subnuclear densities and low temperatures where 
nuclei are so closely packed in a gas consisting mainly of 
dripped neutrons as to melt into uniform matter, nuclear 
matter is expected to possess spatially periodic struc- 
tures composed of rodlike and slablike nuclei and rodlike 
and spherical bubbles 0; these nuclei and bubbles are 
often referred to as nuclear "pasta." Seminal works by 
Ravenhall et al. and Hashimoto et al. Q based on 
liquid-drop models show that a subtle balance between 
nuclear surface and Coulomb energies determines the en- 
ergetically most favorable "pasta" shape. Such "pasta" - 
like structures may be encountered in the outer part of 
a stellar remnant of collapse-driven supernova explosion. 
In this part, the inhomogeneous nuclear matter would be 
neutralized and roughly /3-equilibrated by a relativistic 
degenerate gas of electrons and, during and just after col- 
lapse, of electron neutrinos. [Hereafter, such material will 
be denoted as supernova matter or neutron star matter, 
according to whether the degenerate neutrino gas coex- 
ists or not.] The charge screening action of the electron 
gas on clumps of protons reduces the Coulomb energy 
induced by the protons while increasing the kinetic en- 
ergy of the electron gas. These energy corrections, which 
have yet to be examined in detail, may affect the density 
region in which the "pasta" phases are preferred over the 
low density phase with roughly spherical nuclei and the 
high density phase of uniform matter. 

In this paper we ask the question of how the electron 
screening changes the energies of zero-temperature su- 
pernova matter and neutron star matter at subnuclear 
densities. For the purpose of answering this question, it 
is instructive to consider the two opposite limits: per- 
fect screening and no screening. In the case of perfect 



screening where the electrostatic energy vanishes and the 
kinetic energy of the electron gas is maximal, a coexist- 
ing state of a nucleon liquid of density comparable to 
the normal nuclear density no ~ 0.16 fm~ 3 and a low 
density neutron gas can occur at a unique pressure as a 
phase separated state. On the other hand, in the absence 
of charge screening, i.e., in the homogeneous limit of the 
electron gas, the electrostatic energy is maximal while the 
electron kinetic energy is minimal. In the liquid region in 
/? equilibrium, the decrease in the electron Fermi energy 
as compared to the perfect screening case plays a role 
in increasing the proton fraction and thus lowering the 
bulk energy of the nucleon liquid. Due to such decrease 
in the electron and nucleon bulk energy, a liquid-gas co- 
existing state can occur for a finite range of pressures 
as a mixed state composed of alternating liquid-gas re- 
gions. At sufficiently low densities, this mixed state is 
expected to manifest itself as a state in which the liquid 
regions correspond to roughly spherical nuclei and the 
gas regions are almost vacant. At densities just below 
no, as suggested by earlier studies 0, such a mixed state 
may appear as the "pasta" phases mentioned above, de- 
pending on uncertain quantities such as the neutron and 
proton chemical potentials of a neutron gas and the sur- 
face tension in the presence of dripped neutrons. A real 
situation - partial screening - lies between those two lim- 
iting cases, whereas it remains to be clarified how such 
screening affects the density region of the "pasta" phases 
through modifications of the bulk and electrostatic ener- 
gies. 



The key quantity of the electron screening is the ratio 
of the scale of a nucleus or bubble to the Thomas-Fermi 



2 



screening length A^p as given by 
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where ni '' and /li '' are the averaged electron number 

density and chemical potential, respectively. If the ratio 

is much smaller than unity, the electron density can be 

assumed to be everywhere constant to a good approxima- 
te) 

tion. The typical values of A^p at subnuclear densities 
for neutron star matter (upper) and supernova matter 
(lower) can be estimated in the massless limit as 
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where a is the fine structure constant, ki°^ = (37r 2 ni°-') 1 / 3 
is the electron Fermi wave number, and n and x are the 
nucleon density and proton fraction in the nucleon liquid. 

fe) 

These values of A-pp are larger than the typical values of 
half thickness of a nucleus or bubble, ~ 5 fm, and even 
half the internuclear spacing, ~ 10-15 fm (see, e.g., Refs. 
PUSES IE El) Consequently, the standard approxima- 
tion in which no screening is included seems fairly valid. 
As we shall show within the linearized Thomas-Fermi ap- 
proximation [13, the screening correction to the density 
region in which the phases with bubbles and nonspher- 
ical nuclei appear is indeed small for the typical values 
of the surface tension, while the screening makes such a 
density region larger rather than narrower. This feature 
will be concluded to be general by investigating how the 
electron screening affects the conditions for instabilities 
against fission of spherical nuclei and proton clustering 
in uniform matter. 

Quantum molecular dynamics simulations that have 
recently been performed by one of the authors success- 
fully reproduce the "pasta" phases expected to occur 
as the ground state, without imposing any assumptions 
on nuclear shapes |ll|. In this framework, the elec- 
tron screening has not yet been taken into consideration, 
whereas proton screening, i.e., Coulombic polarization of 
the nuclear interior, is automatically incorporated. Be- 
fore examining the influence of the electron screening on 
the simulations, it is useful to focus on the qualitative na- 
ture of the electron screening effects, which will be clari- 
fied in the present work. We remark in passing that the 
proton screening makes the liquid portion deviate only 
slightly from uniformity, although the proton screening 
length is comparable to the scale of the liquid region. 

To investigate the electron screening effects on the 
phase diagram for the ground-state neutron star matter 
and supernova matter at subnuclear densities, we extend 
a compressible liquid-drop model of inhomogeneous nu- 
clear matter immersed in the lepton gas as constructed 



in Refs. Q and Q by incorporating the electron density 
deviation from the non-screening case in the Thomas- 
Fermi approximation up to linear order in the electro- 
static potential [l(J ■ This model, a generalization of the 
one developed by Baym, Bethe, and Pethick [12| (here- 
after BBP), allows us to obtain the phase diagram in a 
way dependent on uncertain values of the surface tension. 
It is essential to take into account such uncertainties; the 
surface tension is related to the equilibrium size of a nu- 
cleus or bubble, which in turn determines the efficiency 
of the screening. 

This paper is composed as follows. In Sec. |n] we con- 
struct a compressible liquid-drop model for nuclei and 
bubbles, and write down equilibrium conditions for zero- 
temperature neutron star matter and supernova mat- 
ter. Expressions for the electrostatic energy including 
the electron screening effects are derived in Sec. Hill The 
influence of the electron screening on the phase diagram 
is discussed in Sec. II VI 



II. MODEL OF MATTER AT SUBNUCLEAR 
DENSITIES 

In this section, we provide the free energy and equilib- 
rium conditions for zero-temperature neutron star matter 
and supernova matter at densities of order 0.1-lno- The 
zero-temperature approximation is considered to well de- 
scribe the energy of matter in the deepest region of neu- 
tron star crusts since the temperature of the crusts is 
typically < 10 9 K, which is much smaller than typical 
nuclear and electronic excitation energies. For densities 
considered here, on the other hand, matter under stellar 
collapse is well characterized by the lepton fraction Yl, 



Yt, 



(0) , 
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where n v is the number density of a uniform gas of elec- 
tron neutrinos, and rib is the averaged baryon density. 
The lepton fraction typically amounts to 0.3-0.38 due to 
the trapping of electron neutrinos [Hi . IT^j . Hereafter, 
we set Yl = 0.3. Since the degeneracy pressure of the 
neutrino gas is much larger than the thermal pressure 
at temperatures of order 10-50 MeV, to a first approx- 
imation, we shall include the energy of the degenerate 
neutrino gas into the energy of zero-temperature neutron 
star matter, and ignore finite temperature effects. 

Following Refs. and @, we consider five phases 
that consist of spherical nuclei, cylindrical nuclei, pla- 
nar nuclei, cylindrical bubbles, and spherical bubbles, re- 
spectively. Each phase is assumed to be composed of a 
Coulomb lattice of a single species of nucleus or bubble 
at a given baryon density n b . We adopt the Wigner-Seitz 
approximation in evaluating the lattice energy with suf- 
ficient accuracy 0- In this approximation, a cell in the 
bec lattice, including a spherical nucleus or bubble of ra- 
dius rN, is replaced by a Wigner-Seitz cell defined as a 
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sphere having radius r c and the same center. A cylindri- 
cal nucleus or bubble having an infinitely long axis and 
a circular section of radius is taken to be contained 
in a cylindrical Wigner-Seitz cell having the same axis 
and a circular section of radius r c in place of a cell in the 
two-dimensional triangular lattice. For a planar nucleus 
with thickness 2rN, a Wigner-Seitz cell is identical with a 
cell in the one-dimensional layered lattice, having width 
2r c . The values of r c for these phases are chosen so that 
each Wigner-Seitz cell may have zero net charge. 

Hereafter we shall concentrate on the charge screening 
of proton clumps by the electron gas. We shall ignore 
curvature effects, nucleon pairing effects, shell effects in 
inhomogeneous nuclear matter , and fluctuation- 

induced displacements of nuclei and bubbles 0>O|, which 
may have significant consequence to the spatial structure 
of nuclear matter. 



A. Energy of matter 

A compressible liquid-drop model for nuclei and bub- 
bles is useful for investigating the nature of the elec- 
tron screening. While we assume neutrons and protons 
to be distributed uniformly inside and outside the nuclei 
or bubbles and electron neutrinos to be an ideal uni- 
form gas as in earlier investigations, we newly include 
the screening-induced deviation, Sn e (r), of the electron 
number density from the unperturbed constant value 
'; 6n e (r) will be explicitly calculated within the linear 
Thomas-Fermi approximation in the next section. Cal- 
culating the contributions up to second order in 5n e (r) 
accurately, we may write the total energy density, Etot, 
averaged over a single cell as 

wn + wl + (1 - u)E n (n n ) + E e [n e (r)] 
+E u (n u ) (nuclei) , 



w N +w L + uE n (n n ) + E e [n e (r)] 
+E u (n v ) 



(bubbles) 
(4) 

Here wn is the energy of the nuclear matter region (the 
region containing protons) in a cell as divided by the 
cell volume, wl is the lattice energy per unit volume, 
n n is the number density of dripped neutrons outside 
the nuclei or inside the bubbles, E n , E e , and E v are the 
energy densities of the neutron matter, of the electron 
gas, and of the neutrino gas, respectively, and u is the 
volume fraction occupied by the nuclei or bubbles: 



n - n n 
n — nb 



(nuclei) , 
(bubbles) 
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where d is the dimensionality defined as d = 1 for slabs, 
d = 2 for cylinders, and d = 3 for spheres, and n is the 
nucleon number density inside the nuclear matter region. 
Note that the contribution of the density deviation 6n e (r) 



is included in Eq. @ through the lattice energy u>l and 
the electron energy E e . 

The expressions for E n , E e , and E v in Eq. Q are given 

by 

E n {n n ) = [W(k n ,0) + m n c 2 } n n , 
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where W(k n ,0) with k n = (37r 2 n n /2) 1/ ' 3 is the energy 
per neutron for uniform neutron matter, V c is the cell 
volume, and k v = (6n 2 n u ) 1 ^ is the neutrino Fermi wave 
number. 

For the energy of the nuclear matter region, u>n, we 
adopt a generalized version of the compressible liquid- 
drop model developed by BBP which gives rise to 
i 

u> N (n, x,n„,r N ,r c ,d) 

un[(l — x)m n + xm p ]c 2 + unW{k, x) 

+w suri (n, x, n n , r N , u, d) + w c (n,x,r N ,u, d) 

(nuclei) , 

(1 — u)n[(l — x)m n + xm p ]c 2 

+ (1 - u)nW(k, x) + w sm f(n, x, n n , r N , u, d) 
+w c {n,x 1 r Nl u,d) 

(bubbles) , 

where m n (m p ) is the neutron (proton) rest mass, 
W(k, x) is the energy per nucleon for uniform nuclear 
matter of nucleon Fermi wave number k — (37r 2 n/2) 1 ^ 3 
and proton fraction x, as given by BBP [see Eq. (3.19) 
in Ref. |12(], u> sur f is the nuclear surface energy per unit 
volume, and wq is the self Coulomb energy (per unit 
volume) of protons contained in a cell. 

This expression for wn includes three parameters Ci, 
C2, and C3, which are associated with uncertainties in 
the proton chemical potential /ip ^ in pure neutron matter 
as contained in W(k, x) and those in the nuclear surface 
tension 



(12) 



The parameter C\ determines the magnitude of fi p ^ (not 
including the rest mass) as [Eq. (4) in Ref. @] 



u (0) 
H'p 



(13) 
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Hereafter we shall set C\ = 400 MeV frn 2 ; this case is 
consistent with the n n dependence of /Zp " 1 obtained from 
various model calculations as exhibited in Fig. 1 of Rcf. 
§. The other two parameters Ca and C3 are defined as 
[Eq. (6) in Ref. 0] 



as well as under charge neutrality, 



xnu (nuclei) , 

xn(l — u) (bubbles) 



(16) 



E sm { = C* 2 tanh ( -Jr 



pBBP 



(14) 



where = dE n /dn n — m n c 2 is the neutron chemical 
potential in the neutron gas (not including the rest mass) , 
and Ef^f is the BBP-type surface tension [Eq. (7) in 
Ref. |f]. As shown in Fig. 1 of Ref. @ and Fig. 1 of Ref. 
0, the surface tension E SUI f calculated for C2 = 1.0 and 
C3 = 3.5 MeV agrees well with the Hartree-Fock results 
obtained by Ravenhall, Bennett, and Pethick using 
a Skyrme interaction. In the present work, we thus as- 
sume C3 = 3.5 MeV, while we give Ci a range of values 
including unity, 0.1-10. To set such various values for 
the parameter C2 determining the strength of -E S urf al- 
lows us to investigate the electron screening effects more 
clearly since the typical size of the nuclear matter region 
approaches the screening length as C2 becomes larger. 



B. Equilibrium conditions 

Zero-temperature neutron star matter with nuclei or 
bubbles of given shape, in its equilibrium, fulfills the con- 
ditions for stability of the nuclear matter region against 
change in the size, neutron drip, 5-decay, and pressur- 
ization (see Section 2.2 in Ref. 8]). These conditions 
arise from minimization of the energy density Etot with 
respect to five variables n, x, n n , nj, and u at vanishing 
n v and fixed baryon density rib given by 



un + (1 — u)n n 
(1 — u)n + un n 



(nuclei) , 
(bubbles) , 



(15) 



In order to obtain the equilibrium conditions for super- 
nova matter, one has only to repeat the minimization 
at fixed Yl, Eq. 1(5} , and rib, Eq. i|15|) . rather than at 
vanishing n v and fixed rib- 

The expression for the size equilibrium that can be 
obtained from optimization of E to t with respect to tn at 
fixed n, x, n n , and it is 



dr N 



(w S urf + wc+l + E e [n e (r)]) 



= 



(17) 



n,x,n n ,u 



This expression does not agree with the non-screening 
formula [see Eq. (14) in Ref. 0], 



w SUTi = 2w c +l , (18) 
since the total electrostatic energy density denoted by 
wc+h = w c + w l is uo longer proportional to for 
fixed u (see Sec. MIL 

The /3-equilibrium condition can be written as 



Me = /4 N) - M P N) + + (m» - m p )c 2 (19) 



where /i p N ^ is the proton chemical potential in the nuclear 
matter region given by Eq. (20) in Ref. fi v = hk v c is 
the neutrino chemical potential, which vanishes for neu- 
tron star matter, and 
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9 n(l — u)V c dx 



5n e (r) 
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(nuclei) , 
(bubbles) 



(20) 



is the electron chemical potential. 

The expressions for the remaining two equilibrium con- 
ditions are unchanged from the non-screening case 0, Q . 
The condition for the drip equilibrium reads 



(21) 



where /il N ^ and are the neutron chemical potentials 



in the nuclear matter region and in the neutron gas given 
by Eqs. (16) and (17) in Ref. @, respectively. The pres- 
sure equilibrium condition can be expressed as 



p(N) = p(G) 



(22) 



where P< N) and P( G > are the pressures of the nuclear 
matter region and of the neutron gas given by Eqs. (22) 
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and (23) in Ref. @, respectively. 

In order to see the difference from the non-screening 
case explicitly, it is instructive to write down the expres- 
sions for and f»( N > that can be obtained after ex- 
pressing x, rj$, and r c in terms of rib, n, and n n through 



the baryon density i|15|) . the size equilibrium condition 
(|17|l . and the /3-equilibrium condition 119|l . The neutron 

chemical potential p,^ in the nuclear matter region reads 



/' 



(N) 



W(k,x) + 



W(k,x) + 



k dW(k, x) 1 dwc+i 
3 dk + 



d dE sul -f 



u dn 

kdW(k,x) 1 dwc+L 



tn dn 

d u dE sm{ 



+ x\p e - fjb v - (m n - m p )c 2 



(nuclei) , 



3 dk 1 — u dn 

The pressure of the nuclear matter region reads 



+ 



"\p e — [iv — (m n — m p )c 2 ] (bubbles) . 

(23) 



P 



(N) _ 



f nk dW(k,x) dn dE surf 

„ E su[ { 4- 

3 a/e tn 



dk 



E, 



dn u dE sur f 



surf 



tn 1 — u dn 



rN__d_ 

du dr^i 

x,n n 

r 



(w c +l + E e [n e (r)}) 



ru_d_ 

du drN 



{wq+l + E e [n e (r)]) 



(nuclei) , 
(bubbles) 



h(1 — ii),x,r c 



(24) 



For neutron star matter, these expressions are modified 
from the non-screening expressions, Eqs. (24) and (26) 
in Ref. |8(, due to the difference in the size equilibrium 
between Eqs. (JTZJ and (TSJ|. 

Finally, we can calculate the equilibrium values of 
n, x, n n , tn, and r c and thus the optimal value of E to t 
for neutron star matter (supernova matter) at given rib 
and nuclear shape (at Yl = 0.3 and at given rib and nu- 
clear shape) by incorporating Eqs. I|23|) and (|24|) into the 
conditions for drip and pressure equilibria, Eqs. 121[) and 
(|22|l . By comparing the optimal values of E to t obtained 
for the five crystalline phases and uniform matter, we 
can determine the phase giving the smallest energy den- 
sity at various values of rib and Ci and thereby draw the 
phase diagram for the ground state neutron star matter 
and supernova matter, as will be shown in Sec. IIVI 

III. ELECTROSTATIC ENERGY 

We proceed to derive analytic formulas for the electro- 
static energy density u>c+l in the five crystalline phases 
by taking into account the electron screening within the 
framework of the linearized Thomas-Fermi approxima- 
tion [lfj. In this framework, the energy density func- 
tional e[rii(r), 0(r)] (i — n,p, e, v\ (f>, the electrostatic po- 
tential), which is related to the average energy density 
E tot as E tot = y c _1 / coll d 3 re[n i (r),0(r)], is expanded up 
to second order in i5n e (r). Minimization of the expanded 
energy density functional with respect to 4>(r) and Sn e (r) 
leads to the Poisson equation, 

V 2 0(r) - k 2 J(t) 



= -47rn Q (r) 

_ \ -47re[nx6»(r N - |r|) - ni 0) ] (nuclei) , , , 
~ \ -4ne[nx6(\r\ - r N ) - ni° } ] (bubbles) ' 

where riQ is the local charge density in the non-screening 
limit, and the relation between Sn e (r) and <f>, 

8n e (r) = £^ecb(r) . (26) 

For planar configuration, we take Cartesian coordinates 
in which the z-axis is normal to the nuclear surface and 
the origin is located on the central plane of the nucleus; 
for cylindrical configuration, cylindrical coordinates in 
which the radial coordinate p is normal to the surface of 
the nucleus or bubble and the line of p — coincides with 
the symmetry axis of the nucleus or bubble; for spherical 
configuration, spherical coordinates in which the radial 
coordinate r is normal to the surface of the nucleus or 
bubble and the origin is located at its center. In physical 
realizations, the electrostatic potential and its derivative 
must be continuous throughout the system. Thus, appro- 
priate boundary conditions are that the derivative of the 
electrostatic potential at the origin and the cell bound- 
ary is zero, i.e., 4>'(0) — and 4>'(r c ) — 0; and that the 
electrostatic potential and its derivative are continuous 
at the surface of the nucleus or bubble. 

We turn to the solutions to the Poisson equation l|25() 
for slablike nuclei, cylindrical nuclei and bubbles, and 
spherical nuclei and bubbles. The results for the electro- 
static potential read 
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• slab 



± 



Aimxe 



± 



Aimxe 



sinh[K e (r c - r N )] 
sinh(K e r c ) 

sinh(« e rN) 



cosh(«; e |z|) + (1 — u) 



cosh[K e (r c — \z\)\ — u 



• cylinder 



±- 



±- 



Aimxe 
Aimxe 



sinh(« e r c ) 



0< \z\ <r N , 
r c > \z\ > r N , 



(27) 



^i(K e r c ) 

L -'l( K e r cJ 



Ke^N-folKeP) + (1 - It) 



< p < r N , 



(28) 



• sphere 



± 



Aimxe 



±- 



47rnxe 



^e(r c - r N ) cosh[K e (r c - r N )] + (K e r c K e r N - 1) sinh[K e (r c - r N )] sinh(K e r) 
K e r c cosh(/-c e r c ) — smh(ft e r c ) n e r 

K e r^ cosh(K e r^) — sinh(K e rN) n e r c cosh[K e (r c — r)] — sinh[K e (r c — r)] 



+ (!-«) 



K e r c cosh(/-c e r c ) — sinh(K e r c ) 



r 



< r < r N , 

r c > r > r N , 
(29) 



Correspondingly, the electric field E can be obtained 



where the upper plus (lower minus) sign corresponds to 
nuclei (bubbles); /„ and K n are the n-th modified Bcsscl from the derivative of the potential <j) as 
functions. 

I 



• slab 



E(z) 



Aimxe sinh[ft e (r c - r N )] z , 

± l —- smh(K e \z ) i-, < z < r N , 

K e smh(K e r c ) \z\ 

Aimxe sinh(K e rN) . , , -2 

± t-tt r sinh K e (r c - \z ) ^-r r c > z > r N , 

K e sinh K e r c z 



(30) 



cylinder 



' Aimxe 7i(«; e r c )i^i(/t e r N ) - /ri(K e r c )7i(/t e r N ) 



± 



£(P) = < 



±- 



Aimxe 



K e r N /i(K e r N ) 



Ii(K e r c ) 

Ki(n e r c ) 
Ki(K e p) — -Ii(n e p) 



K e r N /i(/t e p) 0<p<r N , 
r c > p> r N , 



(31) 



sphere 



£(r) = ^ 



± 



± 



h(K e r c ) 

Aimxe K e (r c - r N ) cosh[K e (r c - r N )] + (K e r c n e r N - 1) sinh[K e (r c - r N )] 
K e K e r c cosh(«; e r c ) — sinh(«; e r c ) 

K e rcosh(K e r) — sinh(« e r) 
{n e r) 2 

Aimxe K e r^ cosh(K e r^s) — sinh(K e rN) 
K e K e r c cosh(K e r c ) - sinh(K e r c ) 

K e (r c ~ r) cosh[/t e (r c - r)} + {n e r c n e r - 1) sinh[K e (r c - r)] 

(K e r) 2 



< r < tn, 



(32) 



r c > r > r N , 
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where the upper plus (lower minus) sign corresponds to 
nuclei (bubbles) as in Eqs. (|77jl - TM . 

We remark that the electron density deviation Sn e , 
which can be obtained by substituting the solutions for 
4> given by Eqs. (|23)— (EHJ into Eq. does not violate 
charge neutrality in a cell. This is because J n Sn e d 3 r oc 

^11^=0. 

The total electrostatic energy Wc+l in a cell is given 

by 



W c+ L = 7^ / E'd'r . 



8tt 



The electrostatic energy density wq+l can thus be cal- 
culated by substituting Eqs. (J201 E3) into Eq. (33Jl as 



o I \2 2 /.(screen)/ \ 



(33) with 



(34) 



cell 



slab 



s sinh £ 1 sinh t , . , sinh i 

5— H — — smh(s -t) + 2 , cosh(s - 1) - 1 

£ sinh s t sinh s sinh s 



(35) 



• cylinder 



/ 2 (screen) (M) 



= m 



Kx(t) r Ji(t) 



A(t) Ik 



7 ) ^?(^2 



1 K^s) 1 



'^K Q {t)K 2 {t) - Mll (t)I 2 (t) 2^j- - /o(*)/a(t))} 

J («)I 2 («) ^o(s)^ 2 (s) 



A 2 (s) 



K 2 ( S ) 



i 2 A 



r*- 



i ' 2 

1, 2, 0, 



^22 

"91 



x > 2 

1, 2, 0, 



(36) 



sphere 



/i scrcon) (M) 



2t 6 (s cosh s — sinh s) 2 

x j- (l + (s - t)(s 2 t -s- t)) + (l + s(s - *)) cosh[2(s - t)] 

+ i (s 2 t - 4s + t) sinh [2 (s - <)]| (i coshi - sinht) 2 

+ {(s-t) cosh(s -t) + (st - 1) sinh(s - t)} 2 (t 2 + tcosht sinhf - 2 sinh 2 i 



(37) 



where G 2 | is the Meijer's G-function defined as 

1; | 0, o) = ^ j h{z)K 1 {z)zdz . (38) 
Figs, n and [3 we plot the ratio 



^.22 
U 24 



In 



(screen) 



{hi e r c ,n e r^)/ fd{u), where is the non- 



screening limit given by 



fd(u) = 



1 



d+2 



d-2 



du l - 2 ' d 



(39) 



This expression can be obtained by taking K e — > in Eqs. 
(13511 13 7|) . We find from Fig. ^ that as K e increases (i.e., 
the electron density ni°^ increases), the ratio j^ scrocn ) / f d 
decreases most rapidly for the slablike configuration while 



most slowly for the spherical configuration. This configu- 
ration dependence holds also for increase in and r c , as 
can be seen in Fig.|3 We may thus conclude that for fixed 
values of tn , r c , and , the screening-induced reduc- 
tion in the electrostatic energy is larger for the structure 
of lower dimensionality. 

In Figs. |3] and we illustrate the charge distributions 
for the five crystalline phases in neutron star matter and 
supernova matter, respectively. We can see from these 
figures that in both cases the magnitude of the electron 
density deviation \5n e \ induced by the screening is much 
smaller than that of the unperturbed charge number den- 
sity |n<g|/e. This ensures the validity of the linearized 
Thomas-Fermi approximation adopted here. 
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r c =10 fm 




K e (fm- 1 ) 

FIG. 1: Contour plot of the ratio f d SCICcn ^/fd on the K e versus 
tn plane for (a) slablike nuclei (d = 1), (b) cylindrical nuclei 
or bubbles (d — 2) , and (c) spherical nuclei or bubbles (d — 3) . 
The Wigner-Seitz cell radius r c is fixed at r c = 10 fm. In the 
Coulomb limit of K e — > 0, the ratio reduces to unity. 



IV. SCREENING CORRECTIONS 

We now examine the influence of the electron screening 
on the phase diagrams of the ground-state neutron star 
matter and supernova matter at subnuclear densities. In 
Figs. and H3 we draw the phase diagrams of neutron 
star matter and supernova matter, respectively, on the 
versus C2 plane for the cases with and without screening. 
As can be seen from these figures, the screening leads 
to slight expansion of the density region in which the 
"pasta" phases appear, and this expansion is larger in 




1 2 1 2 



K e In 

FIG. 2: Contour plot of the ratio /i scre / fd on the K e 7"N 
versus K e r c plane for (a) slablike nuclei (d — 1), (b) cylindrical 
nuclei or bubbles (d — 2), and (c) spherical nuclei or bubbles 
(d = 3). 



supernova matter than in neutron star matter. 

In order to consider why the screening is more effec- 
tive in supernova matter than in neutron star matter, it 
is useful to plot the screening length A^p and the spatial 
scales, tn and r c , as in Figs. [7| and 0D By comparing these 
figures, we can see that the ratios ^n/A^f and r c /A^ F , 
whose increase makes the electron screening more effec- 
tive, is larger in supernova matter than in neutron star 

matter. This is because of smaller A^ F in supernova mat- 
ter. 

We also see from Figs.[7|and|S]that the electron screen- 
ing acts to increase tn and r c by a small amount, while 
keeping the ratio rjq/r c , or equivalently, u, almost un- 
changed. Such increase in tn stems from the fact that as 
discussed in Sec. [n] the electron screening modifies the 
condition for size equilibrium from Eq. I|18fl to Eq. I|17|l . 
which can be rewritten as 

Wsurf = 2(w c+Lj0 + Sw c +h + SE e ) . (40) 

Here, u>c+l,o = 27r(nxe) 2 r^,u/ £ ;('u) is the total electro- 
static energy density in the non-screening limit, and 
Swq+l and SE e are the corrections to Eq. due to the 
screening-induced changes in the electrostatic energy and 
the electron energy, respectively, which are proportional 
to Kg in leading order. Up to 0(/tg), the sum Swc+l and 
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p .„ hole : = 
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FIG. 3: Charge distribution in a Wigner-Seitz cell in neu- 
tron star matter for the phases with (a) spherical nuclei, (b) 
cylindrical nuclei, (c) slablike nuclei, (d) cylindrical bubbles, 
and (e) spherical bubbles, calculated at typical baryon den- 
sities for C2 = 1.0. The solid lines denote the unperturbed 
local charge density nq(r), divided by e, and the dashed lines 
denote the screening- induced deviation Sn e (r) of the electron 
number density, multiplied by ten. The vertical dotted lines 
are the cell boundaries. 



FIG. 4: Charge distribution in a Wigner-Seitz cell in su- 
pernova matter for the phases with (a) spherical nuclei, (b) 
cylindrical nuclei, (c) slablike nuclei, (d) cylindrical bubbles, 
and (e) spherical bubbles, calculated at typical baryon den- 
sities for C2 = 1.0 and YL = 0.3. The solid lines denote the 
unperturbed local charge density nQ(r), divided by e, and the 
dashed lines denote the screening-induced deviation Sn e (r) of 
the electron number density. The vertical dotted lines are the 
cell boundaries. 



SE e is negative since Swc+l, which is negative, is more 
important than 5E e , which is positive. Consequently, the 
screening makes the equilibrium value of tn larger than 
that in the non-screening limit, given by 



„(°) 



dE. 



surf 



4ir(nxe) 2 f d 



1/3 



(41) 



On the other hand, the negligible screening effect on u 
suggests that the pressure corrections due to the screen- 
ing through the Coulomb pressure and the electron pres- 
sure [see the last two terms in the right side of Eq. Q24[l] 
are negligibly small. 

We proceed to see how the phase structure changes 
with the strength of the surface tension, C*2, and the 
baryon density rib- We find from Figs. [5] and that 
as C2 decreases, the phase boundaries in the case with 
screening approach those in the case without screening. 
This is consistent with the fact that for weaker surface 
tension, the equilibrium size of the spatial structure be- 
comes smaller, leading to smaller r^j/A^p and r c /A^ F . 
As rib increases with C2 fixed, on the other hand, the 
screening induced change in the phase boundaries be- 
comes more appreciable; the increase in the transition 
density between the cylindrical hole to the spherical hole 
phase is larger than that between the slab to the cylindri- 
cal hole phase while being smaller than that between the 
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\ \ hole 
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3 0.1 0.12 
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0.14 0. 



FIG. 5: Zero-temperature phase diagram of neutron star 
matter on the rib versus C'2 plane. The solid lines are the 
phase boundaries obtained for the case allowing for the elec- 
tron screening. The dashed lines are for the case ignoring the 
electron screening, which are taken from the lower left panel 
in Fig. 3 in Ref. ||]. 
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spherical 
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FIG. 6: Zero-temperature phase diagram of supernova mat- 
ter on the rib versus C2 plane. The solid lines are the phase 
boundaries obtained for the case allowing for the electron 
screening. The dashed lines are for the case ignoring the elec- 
tron screening,which are taken from the lower left panel in 
Fig. 3 in Ref. 0. 



spherical hole to the uniform phase. This is partly be- 
cause A^p decreases with increasing density and partly 
because at fixed rib, the higher dimensionality has the 
larger equilibrium values of tn and r c . This dimension- 
ality dependence, which was also obtained in earlier in- 
vestigations based on various nuclear models (see, e.g., 
Refs. 0, IE IE 0)' stems from the fact that generally 
the equilibrium values of the surface energy density w SUI f 
and surface tension E sm { are almost degenerate among 
the five crystalline phases at fixed rib and thus the equi- 
librium value of behaves roughly as r-^ oc d [see Eq. 

COM. 

In order to examine the influence of the screening on 
the phase boundaries in further detail, we list in Table H| 
the transition densities calculated for neutron star matter 
and supernova matter at a typical value of C2 = 1.0. We 
thus find that the "pasta" phases, as a whole, are slightly 
enlarged by the electron screening in such a way that 
the phase boundaries associated with the bubble phases 
move into a higher density, whereas the others move into 
a lower density. Such movement stems from the dimen- 



sionality dependence of the quantity, / 



(screen) 



/ f d , charac- 



terizing the screening correction to the electrostatic en- 
At fixed tn, r c , and ni° , as was found from Figs. 

/fd is smaller and hence the screening 



ergy: 
□ and [2 f 



(screen) 



is more efficient for lower dimensionality. We remark in 
passing that the screening correction to the electrostatic 
energy dominates over the screening correction to the 



- with screening - 
- no screening 




0.08 0.1 

n b (fm- 3 ) 

FIG. 7: Size of a nucleus or bubble, rpj, and of a Wigner- 
Seitz cell, r c , in neutron star matter calculated for C2 = 0.1, 
1.0, and 2.5. The Thomas-Fermi screening length A^p is also 
plotted. The solid lines are the results for the case with screen- 
ing, and the dashed lines are the results for the case without 
screening, which are taken from Fig. 5 in Ref. 8]. The sym- 
bols SP, C, S, CH, and SH stand for sphere, cylinder, slab, 
cylindrical hole, and spherical hole, respectively. 



TABLE I: Transition densities (in fm ) in neutron star 
matter (NSM) and supernova matter (SNM), calculated at 
C2 = 1.0. The lepton fraction YL is set at 0.3 for supernova 
matter. The symbols O an( A x denote the cases with and 
without screening, respectively. The symbols SP, C, S, CH, 
SH, and U stand for sphere, cylinder, slab, cylindrical hole, 
spherical hole, and uniform matter, respectively. 





screening 


SP^C 


C^S 


S^CH 


CH^SH 


SH^U 


NSM 


X 


0.07854 


0.09537 


0.11636 


0.12191 


0.12571 


NSM 





0.07849 


0.09535 


0.11648 


0.12206 


0.12592 


SNM 


X 


0.03433 


0.05364 


0.08965 


0.10364 


0.11812 


SNM 





0.03407 


0.05334 


0.08990 


0.10410 


0.11974 



electron energy. 

The decrease in the transition density between the 
phases with spherical and with cylindrical nuclei can be 
understood from the condition for fission instability of 
spherical nuclei in a Coulomb lattice. Up to C^w 1 / 3 ), 
this condition reads 



uuc > 2w, 



surf 



(42) 



where wq = ATr(nxe) 2 ur^/5 is the Coulomb self energy 
of the nucleus, divided by the cell volume. Condition 
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FIG. 8: Size of a nucleus or bubble, rpj, and of a Wigner- 
Seitz cell, r c , in supernova matter calculated for Yl. = 0.3 and 
C2 = 0.1, 1.0, and 2.5. The Thomas- Fermi screening length 
A^ F is also plotted. The solid lines are the results for the case 
with screening, and the dashed lines are the results for the 
case without screening, which are taken from Fig. 5 in Ref. 
The symbols SP, C, S, CH, and SH stand for sphere, cylinder, 
slab, cylindrical hole, and spherical hole, respectively. 



(|42|l provides a critical size, 



r N 



30E, 



surf 



47r(nxe) 5 



1/3 



(43) 



In the absence of screening, the equilibrium nuclear size 
given by Eq. (|41() in which fd=3 is retained up to 
0(u 1 /' i ), reaches this critical size when u = 1/8. In the 
presence of screening, as in Figs. and 00 the equilibrium 
value of tn is larger than r^' at given n^. Since the 
equilibrium value of tn increases with nb for spherical 
nuclei, the density region in which the spherical nuclei 
are unstable against fission extends to a lower density. 

The increase in the density at which the system turns 
into uniform matter is the most salient feature induced 
by the screening; as can be seen from Table [U the density 
increase for typical values of the surface tension (C2 — 1) 
amounts to ~ 2 x 10 -4 fm~ 3 for neutron star matter 
and ~ 0.0015 fm~ 3 for supernova matter. This feature 
can be understood from the condition for proton cluster- 
ing instability in /3-equilibrated uniform nuclear matter, 
which was originally obtained by BBP |T3] for neutron 
star matter. (For supernova matter, one can follow the 
same argument as that for neutron star matter since the 
presence of trapped neutrinos makes no difference in the 
argument.) They expanded the energy density functional 
E[rii(r)] (i = n, p, e) in the extended Thomas- Fermi 



model with respect to small proton density fluctuations 
6n p (r) around the homogeneous state. The variation of 
the total energy caused by the inhomogeneity yields 



E - En = I 



d 3 q 

(2^)3 



v(q)\6n p (q)\' 



(44) 



where En is the energy of the unperturbed state, <5n p (q) is 
the Fourier transform of Sn p (r), and v(q) is the effective 
potential between protons given by 



v (q) = Vo + I3q 2 + 



Aire 2 
q 2 + k 2 



(45) 



The first term in the right side of Eq. I|45|l is the bulk 
contribution, and the second is the gradient contribution 
(for detailed expressions for vq and (3, see Ref. The 
potential v(q) takes a minimum value t> m i n at q = Q, 
where 



Q 



Aire 2 



1/2 



K e 1 



w + 2(47re 2 /3) 1/2 -/3K c 



(46) 
(47) 



The condition for the proton clustering instability, v m { n < 
0, is satisfied at densities below a critical density ni nst . 
This is because the bulk contribution vq, which is domi- 
nant in i> m i n , is an increasing function of at densities 
around v , f m i n ~ (see, e.g., Fig. 2 in Ref. 0). In 
the presence of the electron screening the critical den- 
sity Hinst becomes higher since the Coulomb interaction 
between charge inhomogeneities, which tends to suppress 
the proton clustering, is weakened by the electron screen- 
ing through the term —f3n e 2 (< 0) in Eq. (|4"7|l . 

In summary we have examined the electron screening 
in the inhomogeneous phases of nuclear matter and clar- 
ified its influence on the zero-temperature phase diagram 
of neutron star matter and supernova matter at subnu- 
clear densities. We have found that the density region 
occupied by the "pasta" phases is slightly expanded by 
the electron screening in such a way that the three phase 
boundaries associated with the bubble phases move into 
a higher density, whereas the other two move into a lower 
density. This expansion stems from the model indepen- 
dent features that the electron screening makes spherical 
nuclei more subject to fission by enlarging the equilib- 
rium size through modifications of the size equilibrium 
condition as in Eq. I|40|l and that the role played by 
the Coulomb interaction between small charge inhomo- 
geneities in uniform nuclear matter in suppressing the 
proton clustering instability is weakened by the electron 
screening. 

Finally, we consider on what physical conditions a 
mixed state composed of two phases having different 
charge density appear as in the nuclear liquid-gas mixed 
state studied here. As mentioned in Ref. [ljj in the 
context of a quark-hadron mixed state, a mixed state 
can be stable only when the charge screening and sur- 
face tension are sufficiently weak. In the present nuclear 
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case, the charge screening is weak and mainly reduces the 
electrostatic energy, leading to a larger density range of 
the mixed state for typical values of the surface tension. 
In the case of the quark-hadron mixed state, the charge 
screening can be strong enough to raise the bulk energy 
of the system and hence prefer the phase separated state 
over the mixed state for realistic values of the surface 
tension, as discussed in Ref. |20|. 
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